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A closed form of the electrostatic potential of a homogeneously charged cube is derived by 
integration. The exact result is compared with multipole expansions for the exterior and interior 
of the cube. The electrostatic potential of a homogeneously charged square in two-dimensional 
electrostatics is also determined. 



I. INTRODUCTION 



The electrostatic potential of a homogeneously charged cube appears in theoretical studies of Wigner lattices [|J . 
In computer simulations of ionic systems using minimum-image electrostatics, it determines the electrostatic self- 
interaction of ions [^| H . In Rcf . Q , Hummer et al. presented a simple calculation of the electrostatic potential at the 
center of a homogeneously charged cube. In this work, a closed form of the electrostatic potential will be determined 
for arbitrary positions. This analytic form can be used for the evaluation of lattice sums. It can also be applied 
as a correction when electrostatic potentials are calculated on a grid, assuming that the grid volumes are uniformly 
charged rather than carrying a point charge at the center. The analytic form of the potential will be compared with 
■ multipole expansions IfO 
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II. CALCULATION OF THE ELECTROSTATIC POTENTIAL OF A CUBE 



The electrostatic potential 4> c of a cube [—1/2, 1/2] 3 with charge density one will be calculated by integration. The 
^ en potential at a point with Cartesian coordinates [u, v, w) can be written as 
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4> c (u, v,w) — I dx j dy I dz [(x — u) 2 + (y — v) 2 + (z — w) 2 ] ^ 2 , (1) 
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" where Gaussian units are used. 4> c can be rewritten as 
d ' 

1/2-it) 

dz(x 2 +y 2 +z 2 y 1/2 . (2) 

l/2-w 

Summation of the results of partial integration with respect to x, y, and z yields a reduction to three two-dimensional 
integrals, 
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2(f> c (u, V, w) = / dx dy z {x 2 + y 2 + z 2 ) 

J-l/2-u J-l/2-v L ' Jz=-l/2-u> 

+ cyclic permutations (x, u; y, v; z, w) — > (y, v; z, w; x, u) and (z, w; x, u; y, v) . (3) 
The two-dimensional integrals can be further reduced using 
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where xq, x\, yo, and y\ are arbitrary integral boundaries. The remaining one-dimensional integrals can be calculated 
using partial integration and conventional substitution for algebraic integrands, 
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dxln[(x 2 + a 2 ) 1/2 + b] = x ln[(x 2 + a 2 ) 1/2 + b] - x 
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+ (x 2 + a 2 ) 1 / 2 + b 
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where 



. , , J arctan(x) for a 2 >b 2 
1 artanh(a;) for a 2 < b 2 . 

Combining the previous results yields a closed form for the electrostatic potential of a unit cube: 
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3 (u, v,w) = - \- c 4 ,i c u+1 
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The integration boundaries are defined as co,o = — 1/2 — u, 01,0 = 1/2 — w, co,i = —1/2 — v, ci,i = 1/2 — w, 
co,2 = —1/2 — ui, and ci,2 = 1/2 — w. The values of I + 1 and Z + 2 in Eq. (0) are defined modulo 3, i.e., co,3 = co,o 
etc. The arctan function to be used in Eq. (0) takes into account the sign of numerator and denominator and yields 
results between -vr and tt ("atan2" in FORTRAN and C). 

An immediate consequence of Eq. (^) is the electrostatic potential at the center of a unit cube 



</) c (0,0,0) = 31n(3 1 / 2 + 2) -| . 
Previous calculations of </> c (0, 0,0) involved rather elaborate manipulations 
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III. CALCULATION OF THE ELECTROSTATIC POTENTIAL OF A SQUARE 



In two-dimensional electrostatics, the charge interaction (Green's function of the Laplacian) is given by — lnr, 
where r is the distance. The electrostatic potential <fi s of a square [—1/2, 1/2] 2 with unit charge density will again 
be calculated by integration. cf> s is also the electrostatic potential of a square cylinder that is infinitely extended in z 
direction. The potential at a point with Cartesian coordinates (it, v) is written as 

I [1/2 rl/2 

<t> s (u,v) = -- dx dy\n[{x - u) 2 + {y - v) 2 } . (9) 

z J -1/2 J -1/2 

Elementary integration yields 

1 1 1 
<t> s {u,v) = -t;J2Y1 

i=0 j=0 

(10) 

where xo = — 1/2 — u, x\ = 1/2 — w, yo = —1/2 — v, and y\ = 1/2 — v. The appropriate arctan function to be used in 
Eq. @, yields values between -tt/2 and tt/2 ("atan" in FORTRAN and C). 



Xi yj in [x^ + yj) — o Xi yj + yj arctan hJj arctan — 
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IV. MULTIPOLE EXPANSION 



The electrostatic potential of a cube can be expanded in "kubic" harmonics, i.e., harmonic functions with cubic 
symmetry ]l|P,p|-pO| . For the exterior, one obtains 

M*) = ~ + Ci K 4 (r) r- 9 + C 6 K 6 (r) r~ 13 + • • • , (11) 
r 

where r = (u, v, w), r = |r|. With T n = u n + v n + w n , the kubic harmonics of order 4 and 6 can be written as (l(J 

K A {v) = T 4 - | r 4 (12) 

if 6 (r)=T 6 -^T 4 r 2 + ^r 6 . (13) 



For this form, the expansion coefficients are C4 = —7/192 and Cq = 11/192 ]l|,|7j. 

For the interior of the cube, we derive the multipole-expansion coefficients of order 2, 4, and 6 from a direct 
Taylor expansion in x direction. The angular dependence can then be inferred by cubic symmetry^] The electrostatic 
potential on the x axis can be expressed as 

/1/2-u 
dx f(x) , (14) 
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Taylor expansion of 4> c {u, 0, 0) around u = yields the expansion 
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Figure ^ shows the electrostatic potential cf> c along the directions (u, 0,0), (u,u, 0), and (u,u,u) calculated from the 
exact result Eq. ^ and the expansions Eq. ([ll]) and ([l6|), both including terms up to K4. The expansions show 
the largest disagreement near the surface of the cube (u = 1/2) where they start to diverge. Otherwise, they closely 



reproduce the exact potential [Eq. ( 11 ) for r — ► 00 and Eq. ( |16| ) for r — ► 0] . 

The divergent behavior reflects an inherent problem of the near- and far-field expansions. By construction, the 
Laplacians of Eqs. (11) and ([ll]) are a delta function at r — and a constant — 47r, respectively, independent of the 



order of the expansions. The former corresponds to a unit point charge and is correct only outside the cube; the latter 
corresponds to a homogeneous charge density and is correct only inside the cube. 

V. CONCLUSION 

Nijboer and Ruijgrok Q analyzed the difference between the energy per particle in a Wigner lattice and the energy 
of a point charge in the field of the other charges. These authors studied an infinite replication of neutral cubes 
consisting of a unit point charge at the center and a compensating background. A reduction of the electrostatic 
potential <fic of a homogeneously charged cube to a one-dimensional integral resulted in 



1 Some higher-order kubic harmonics are degenerate ||, requiring two independent expansion directions to get the correct 
angular dependence. 
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..,(«.<•.«") = - / dit- 2 ^[fr(t*,t)fc(v,t)fc(u;,t)] , (17) 



where 

h(x, t) = erf 
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and erf is the error function.^] The solution of the one-dimensional integral in Eq. ( |l7| ) would give the closed form 
Eq. (Q) of this work. Eq. (Jfy, numerical integration of Eq. (|l7|), and direct Monte Carlo integration of Eq. (|l|) were 
compared for a few hundred points and gave identical results within the error margins of the numerical integration 
in Eq. jl7|) and the statistical errors of the Monte Carlo procedure. Eq. (Q) has the advantage of being analytical. It 
can be evaluated fast and with arbitrary precision on the computer. 

The electrostatic potential <^ c (0, 0, 0) at the center of the cube as listed in Eq. (|J) can be used to correct effectively 
for finite-size effects in computer simulations of ionic systems under periodic boundary conditions, when minimum- 
image electrostatics is used [0. An example is the calculation of single- ion chemical potentials where the 
electrostatic energy of an excess ion has to be calculated. The system-size dependence is greatly reduced if the excess 
charge is compensated with a homogeneous background. The electrostatic energy u of the excess charge q at r = is 
then the sum of the interactions with the other charges qi at r, and with the background, 
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u = qJ2^/n + q 2 M0,0,0)/L , (19) 

i=l 

where a cubical box of length L is used. 

Another application is the calculation of electrostatic potentials when charges are given on a grid, for instance, when 
ionic density distributions are known ]1^| . Usually, the grid charges are assumed to be point charges. In an improved 
description, the charges are smeared out over the grid cells. The electrostatic potentials can then be calculated using 
Eq. (Q) or the multipole expansions Eq. (|ll|) and ([l6|). This eliminates the singularities in the electrostatic potential 
and gives a more accurate description near local charge concentrations. 
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FIG. 1. Electrostatic potential <j> c (u,v,w) of a homogeneously charged cube. Panels A, B, and C show <j> c for directions 
(u, 0, 0), (u, u, 0), and (u, u, u) as a function of u, respectively. The insert in panel A shows a unit cube illustrating the directions 
A, B, and C. The solid line is the exact result Eq. ([?]). The dashed and dot-dashed lines are the expansions Eq. ( pd| ) and ( p^ ) 
up to Ki for the exterior and interior, respectively. The vertical dotted line at u — 0.5 indicates the boundary of the cube. 
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Fig. 1 : G. Hummer 



